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Abstract 

Using data collected by the high energy photoproduction experiment FOCUS 
at Fermilab we performed a Dalitz plot analysis of the Cabibbo favored decay 
iir~7r+7r+. This study uses 53653 Dalitz-plot events with a signal fraction 
of ~ 97%, and represents the highest statistics, most complete Dalitz plot analysis 
for this channel. Results are presented and discussed using two different formalisms. 
The first is a simple sum of Breit-Wigner functions with freely fitted masses and 
widths. It is the model traditionally adopted and serves as comparison with the 
already published analyses. The second uses a K-matrix approach for the dominant 
iS-wave, in which the parameters are fixed by first fitting Ktt scattering data and 
continued to threshold by Chiral Perturbation Theory. We show that the Dalitz 
plot distribution for this decay is consistent with the assumption of two body dom- 
inance of the final state interactions and the description of these interactions is in 
agreement with other data on the Ktt final state. 
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1 Introduction 

In recent years, the Dalitz plot technique has been widely applied to the 
heavy-flavor sector. It is acknowledged as a powerful tool with which to study 
charm and beauty decay dynamics, and so test the consistency of the Standard 
Model, investigate CP violation effects and set limits on new physics. To per- 
form precision studies requires an accurate description of the hadron dynamics 
that colors and shapes the final states. Whether in the three-pion [1], or here 
in the Km: final state, two body interactions with scalar quantum numbers 
dominate the Dalitz plot distribution. Consequently broad overlapping reso- 
nances control the dynamics. A coherent sum of Breit-Wigners can, in general, 
provide an adequate description of data. However, the effective Breit-Wigner 
mass and width need not accurately reflect the true positions of any resonance 
poles, particularly for wide states like the k and a. If three-body interactions 
play a limited role, then wc can adopt a paramctrization which enforces the 
two-body unitarity constraint and is consistent with the two-body scattering 
data. This is naturally embodied in the K-matrix formalism. The FOCUS 
collaboration has already performed a pioneering Dalitz plot analysis of the 
D'^ and Df — > 7r"'"7r~7r''" decays [1] implementing the K-matrix formalism for 
the description of the nn S-wsve intermediate states. It led us to conclude 
that the three pion D-decay is well-described using data from tttt scattering 
and that any cr-like object in D-decay is consistent with the a extracted from 
TTTT scattering [2] . The same conclusion was reached by other authors (see for 
instance [3] and [4]). The D'^ — > K~7r'^7r'^ analysis discussed in this letter 
represents the analogous study in the Kir system. 

It is worth noting that the K-matrix formalism [5,6], originating in the context 
of two-body scattering, can be generalized to cover the case of production of 
resonances in more complex reactions [7] , with the assumption that the two- 
body system in the final state is an isolated one and that the two particles do 
not simultaneously interact with the rest of the final state in the production 
process [6] . The vahdity of the assumed quasi two-body nature of the process 
of the K-matrix approach can only be verified by a direct comparison of the 
model predictions with data. In particular, the failure to reproduce the Dalitz 
plot distribution could be an indication of the presence of relevant, neglected 
three-body effects. 

Within the X-matrix approach, we present the first determination in D-decays 
of the two separate isospin contributions, / = 1/2 and / = 3/2, for the 
S'-wave Kn system. The 1 = 1/2 component is the most important, being 
dominated by broad resonances. In comparison, the / = 3/2 contribution is 
small at low Kn masses, but becomes larger at higher masses. This results in 
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a significant interference between tliese components to describe tlie full Dalitz 
plot distribution. We will see that our results indicate close consistency with 
Kit scattering data, and consequently with Watson's theorem predictions for 
two-body Kn interactions in the low Kn mass region, where elastic processes 
dominate. 



2 Signal selection 

The data for this analysis were collected during the 1996-1997 run of the 
photoproduction experiment FOCUS at Fermilab. The detector, designed and 
used to study the interaction of high-energy photons on a segmented BeO 
target, is a large aperture, fixed-target magnetic spectrometer with excellent 
Cerenkov particle identification and vertexing capabilities. Most of the FOCUS 
experiment and analysis techniques have been described previously [8-10]. 

The FOCUS collaboration has already published the analysis of the meson 
lifetime in the K~7r~^7i~^ channel [11]. Candidates for the present Dalitz plot 
analysis have been selected according to an analogous set of cuts. A decay 
vertex is formed from three reconstructed charged tracks. The momentum 
of the D candidate is used to intersect other reconstructed tracks to form a 
production vertex. The confidence level (C.L.) of each vertex is required to 
exceed 1 %. After the vertex finder algorithm, the variable i, which is the 
separation of the primary and secondary vertices, and its associated error 
ae are calculated. We reduce backgrounds by requiring £/ae > 15. The two 
vertices are also required to satisfy isolation conditions. The primary vertex 
isolation cut requires that a track assigned to the decay vertex has a C.L. 
less than 0.1% to be included in the primary vertex. The secondary vertex 
isolation cut requires that all remaining tracks not assigned to the primary 
and secondary vertex have a C.L. smaller than 0.001 % to form a vertex with 
the D candidate daughters. 

The Cerenkov particle identification is based on likelihood ratios between 
the various stable particle hypotheses [9]. The product of all firing probabil- 
ities for all cells within the three Cerenkov cones produces a x^-like variable 
Wi = —2 In(likelihood) where i ranges over the electron, pion, kaon and proton 
hypothesis. We require AWk = W-,, — Wk > 5 for the kaon candidate, and 
APVjr = Wk — W^TT > 5 for both pions in the final state. Kaon and pion con- 
sistency, AW = Wk - Wrain < 3 and AW = W^, - Wram < 3 is also required 
where VFjnin is the minimum Wi. These Cerenkov cuts reduce the residual con- 
tamination of the Df K^K'^TT^, when a kaon is mis-identified as a pion, 
to a negligible level. The final sample invariant mass distribution is shown in 
Fig. [H Yields for signal and background, evaluated within ±2ct from the mass 
peak, i.e, from 1.85 to 1.89 GeV, consist of 52460 ± 245 and 1897 ± 39 events. 
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respectively. Events satisfying the kinematic limit condition populate the fi- 
nal Dalitz plot, shown in the same Fig. [H on which the amplitude analysis is 
performed. 



^3.25 




Figure 1. K^tt^tt^ Dalitz plot and mass distribution: signal and sideband 

regions are indicated in red and blue respectively. Sideband are at ±(6-8) a from 
the peak. 



3 The Dalitz plot analysis 



Previous analyses of the K'ty^'k'^ used an isobar model, where the 

decay amplitude consists of a simple sum of Breit-Wigner functions. In this 
analysis, the masses and widths for scalar states will be determined by fitting 
to this particular decay channel. This model usefully serves as the standard 
for quality of fit, with no presumption as to the correctness of the physics the 
model embodies. The limitation of this approach, as already mentioned, is that 
the parameters of the scalar states are "ad hoc" , established with no reference 
to those found in other processes, in particular scattering data. Moreover, the 
states involved are wide and overlapping and the Breit-Wigner modelling is 
too simplistic. Assuming dominance of two-body interactions the K-matrix 
formalism [5, 6] provides the correct treatment for the S'-wave Kn system. It 
is instructive to compare the K-matrix against this more commonly employed 
isobar fit and we will present and discuss the results for the two approaches. 
A third approach, namely a partial-wave analysis of the Knir system, will be 
considered in a future publication. 
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3.1 The isobar model 



The decay amplitude in this formahsm is written as a coherent sum of am- 
phtudes corresponding to a constant term for the uniform direct three-body 
decay and to different resonant channels: 

M = aoe'^° + J2 «i e'^' B{abc\r), (1) 
j 



where a, b and c label the final-state particles. Coefficient and phase of the 
K*{S92), our reference amphtude, are fixed to 1 and respectively. 

B{abc\r) = B{a,b\r)S{a,c) where B{a,b\r) is the Breit-Wigner function 



and S{a, c) = 1 for a spin-0 resonance, S{a, c) = —2a • c for a spin-1 reso- 
nance and S{a, c) = 2 (3(a ■ c)^ — |ap|cp) for a spin-2 state. The a and c are 
the three-momenta of particles a and c measured in the ab rest frame. The 
momentum-dependent form factors Fjy and represent the strong coupling 
at each decay vertex, and are of the Blatt-Weisskopf form. For each resonance 
of mass Mr and spin J we use a width 




where p is the decay three-momentum in the resonance rest frame and the sub- 
script r denotes the on-shell values. The order of particle label is important for 
defining the phase convention; here the first particle is the opposite-sign one, 
i.e., the kaon. Each decay amplitude is then Bose-symmetrized with respect 
to the exchange of the two identical pions. 



3.2 The K-matrix model 



The model of 5'-wave states requires particular care in order to account for 
non-trivial dynamics generated by the presence of broad and overlapping res- 
onances: a real K-matrix guarantees unitarity for two-body interactions. An 
additional complication in the Ktt system comes from the presence in the 
S-wme of the two isospin states, / = 1/2 and / = 3/2. 
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The K-matrix form we use as input describes the 5'-wave K'-k^ K~7r~^ 
scattering from the LASS experiment [12] for energy above 825 MeV and 
K~iT~ — >■ K~TT~ scattering from Estabrooks et al. [13]. The K-matrix form 
follows the extrapolation down to Ktt threshold for both J = 1/2 and 1 = 3/2 
S-wave components by the dispersive analysis by Biittiker et al. [14], consistent 
with Chiral Perturbation Theory [15]. The complete form is given below in 
Eqs. (1319]) with the parameters listed in Table [T] [16]. 

Although only the / = 1/2 is dominated by resonances, both isospin compo- 
nents are involved in the decay of the meson into K^-k^-k^ . A model for 
the decay amplitudes of the two isospin states can be constructed from the 2 
X 2 K-matrix describing the / = 1/2 S'-wave scattering in {K'k)i and {Kri')2 
(with the subscripts 1 and 2, respectively, labelling these two channels), and 
the single-channel K-matrix describing the 7 = 3/2 K~t[^ K'tt'^ scatter- 
ing. 

The total D decay amplitude in Eq. ([1]) can be written as 

M = {Fy,)^{s)+Fy,{s) + Y.aj e''^ B{abc\r), (4) 

j 

where s = M'^{K-n), (Fi/2)^ and F3/2 represent the J = 1/2 and / = 3/2 
decay amplitudes in the Kn channel, j runs over vector and spin-2 tensor 
resonances 0, and B{ahc\r) are Breit-Wigner forms as in Eq. ([T]) and Eq. ([2]). 
The J > resonances should, in principle, be treated in the same K-matrix 
formalism. However, the contribution from the vector wave comes mainly from 
the i^'*(892) state, which is well separated from the higher mass i^'*(1410) and 
/('*(1680), and the contribution from the spin-2 wave comes from i^'2(1430) 
alone. Their contributions are limited to small percentages, and, as a first 
approximation, they can be reasonably described by a simple sum of Breit- 
Wigners. More precise results would require a better treatment of the over- 
lapping i^"*(1410) and -ft'*(1680) resonances as well. F1/2 is actually a vector 
consisting of two components: the first accounting for the description of the 
Ktt channel, the second of the Krj' channel: in fitting —* K^tt^tt~^ we 
need, of course, the (-Fi/2)^ element. Its form is 

(Fi/2)l = (/-^i^l/2p)l,'(Pl/2)„ (5) 

where / is the identity matrix, K1/2 is the K-matrix for the I = 1/2 S- 
wave scattering in Kn and Kr]', p is the corresponding phase-space matrix 
for the two channels [6] and (Pi/2)j is the production vector in the channel 

^ Higher spin resonances have been tried in the fit with both formalisms but found 
to be statistically insignificant. 
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j. In this model [7], the production process can be viewed as consisting of an 
initial preparation of states, described by the P-vector, which then propagates 
according to (/ — iKp)~^ into the final one. 

The form for F3/2 is 

Fs/2 = {I-tKs/2p)-'P3/2, (6) 

where is the single-channel scalar function describing the 7 = 3/2 

K~7c^ K~7T~^ scattering, and P3/2 is the production function into Ktt. 

The Kj matrix for the isospin / state is derived by fitting scattering data via 
the corresponding 7} matrix defined as 

Tj = {I-iKjp)-'Kj. (7) 

The (^/2)ii and 7^/2 functions are composed to fit the K~7t~^ — ^ K^tt^ S- 
wave through the Clebsch-Gordan coefficients to give Tu = ('7i/2)ii + 73/2/2. 

Fitting of the real and imaginary parts of the K^tt^ K^tt^ LASS ampli- 
tude, shown in Fig. [21 and using the predictions of Chiral Perturbation Theory 
to continue this to threshold, gives the K-matrix parameters in Table [TJ 



LASS LASS 

0.6 1 1 1 1 , 1 1 1 , 1 1 1.2 1 1 1 1 1 1 — 




-0.8 1 1 1 1 ' 1 1 1 ' 1 1 -0.2 1 ' ' 1 1 ' 1 1 ' ' 1 

0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 



M (Kjt) (GeV) M (K;r) (GeV) 

Figure 2. Real and imaginary K^t:^ K^ir^ amplitudes from the LASS experi- 
ment and their K-matrix fit results. 

The 1=1/2 K-matrix is a one-pole, two-channel matrix whose elements are 
given in Eq. ([H]). 
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K22 — ( ( 1" C22O + C221S + C222S^ 



Ku = { ) ( + C120 + Ci2i5 + C122S' ] , (8) 



where the factor of Snorm = i^^'k + conveniently introduced to make 

the individual terms in the above expression dimensionless. gi and g2 are the 
real couplings of the Si pole to the first and the second channel respectively. 
Sgi = 0.23 GeV^ is the Adler zero position in the / = 1/2 ChPT elastic scat- 
tering amplitude [!]• Cm^ C22i and Ci2i for i = 0, 1, 2 are the three coefficients 
of a second order polynomial for the diagonal and off-diagonal elements of 
the symmetric K-matrix. Polynomials are expanded around s = s/snorm — 1- 
This form generates an S-matrix pole, which is conventionally quoted in the 
complex energy plane as E = M — iV /2 = 1.408 — zO.llO GeV. Any more 
distant pole than ^0(1430) is not reliably determined as this simple K-matrix 
expression does not have the required analyticity properties. Nevertheless, it 
is an accurate description for real values of the energy where scattering takes 
place. Numerical values of the terms in Eq. ([8]) are reported in Table [H 

pole (GeV^) coupling (GeV) Cm Ci2i C22i 

si = 1.7919 

gi = 0.31072 
52 = -0.02323 

Clio 
Cm 
C112 

Table 1 

Values of parameters for the 1 = 1/2 K-matrix. 

The 1 = 3/2 K-matrix is given in Eq. ([9]). Its form is derived from a simulta- 
neous fit to LASS data [12] and to K~tt~ — » K~tt~ scattering data [13]. It is 
a non-resonant, one channel scalar function. 



= 0.79299 C120 = 0.15040 C220 = 0.17054 
= -0.15099 C121 = -0.038266 C221 = -0.0219 
= 0.00811 C122 = 0.0022596 C222 = 0.00085655 



^0^ 



^3/2= ( ^1 (/^iio + ^iiiS + 1^1125'). (9) 



^ Chiral symmetry breaking demands an Adler zero in the elastic S-wave ampli- 
tudes in the unphysical region. ChPT at next-to-leading order fixes these positions 
SOI [14,15]. 
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In Eq. ([9]) Sgs = 0.27 GeV^ is the Adler zero position in the 1 = 3/2 
ChPT elastic scattering and the values of the polynomial coefficients are 
Di^Q = -0.22147, Dm = 0.026637, and ^112 = -0.00092057 [16]. 

When moving from scattering processes to /^-decays, the production P-vector 
has to be introduced. While the K-matrix is real, P-vectors are in general 
complex reflecting the fact that the initial coupling — >• {K~ 'K^)'K^pectator 
need not be real. The P-vector has to have the same poles as the K-matrix^ 
so that these cancel in the physical decay amplitude. Their functional forms 
are: 



(Pi/2)i = + (Clo + cns + c^2s')e'^' 



Si — s 



(10) 



(^1/2)2 = + (C20 + C2lS + C22s')e'^' 



:iii 



Si — s 



P3/2 = {c30 + C3iS + C32S^)e'^'. (12) 

/3e*^ is the complex coupling to the pole in the 'initial' production process, 
Qi and g2 are the couplings as given by Table [H The Kir mass squared 
Sc = 2 GeV^ corresponds to the centre of the Dalitz plot. It is convenient to 
choose this as the value of s about which the polynomials of Eqs. fll0til2p are 
expanded, by defining s = s — Sc- The polynomial terms in each channel are 
chosen to have a common phase 7, to limit the number of free parameters in the 
fit and avoid uncontrolled interference among the physical background terms. 
Thus the coefficients of the second order polynomial, qj, are real. Coefficients 
and phases of the P-vectors, except gi and g2, are the only free parameters of 
the ffi determining the scalar components. 



3. 3 The likelihood function and fitting procedure 



In analogy with our previous works we perform a fit to the Dalitz plot 
with an unbinned likelihood function, C, consisting of signal and background 
probability density. The signal probability density is corrected for geometri- 
cal acceptance and reconstruction efficiency. The shape of the background in 
the signal region is fixed to that derived from a ffi to the Dalitz plot of the 
mass sidebands. It is parametrized through an incoherent sum of a polynomial 
function plus resonant Breit-Wigner-like components. The signal fraction is 
estimated by a ffi to the K'Tr'^n'^ mass spectrum. Checks for ffiting proce- 
dure are made using Monte Carlo techniques and all biases are found to be 
small compared to the statistical errors. The systematic errors on our results 
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include split-sample and fit-variant components [1,17]. The sample is split in 
three different ways: low and high D momenta, D'^ and D~ , and according to 
the two main running periods of the experiment. 

The fit-variant systematics is evaluated by varying the background parametriza- 
tion used to fit the sideband Dalitz plot, moving the sideband regions, and 
letting the background freely float within the error returned by the sideband 
fit. 



3.4 Isobar model results 

We allow for the possibility of contributions from all known well-established 
{K~7i^) resonances [18]. In addition, a constant amplitude accounts for the di- 
rect decay of the D meson into a non-resonant three-body final states. The fit 
parameters are amplitude coefficients and phases Si of Eq. ([Tj). Contributions 
are removed if their amplitude coefficients, have less than 2 a significance 
and the fit confidence level increases. The fit confidence levels (CL.) are eval- 
uated with a estimator over the Dalitz plot with bin size adaptively chosen 
to maintain a minimum number of events in each bin. A fit consisting solely of 
well-established K~tt~^ resonances with masses and widths as in the PDG [18] 
results in a very poor solution, with an adaptive bin x^/d.o.f of more than 3, 
and with a very high level of the non- resonant component, about 90% of the 
total, atypically high for charm decays. 

Following a previous work [19] a low mass K~7r~^ resonance k is introduced 
with mass and width freely floating in the fit. The insertion of this state, whose 
mass and width are fitted to be 883 ± 13 MeV/c^ and 355 ± 13 MeV/c^ re- 
spectively, returns a x^/d.o.f of more than 2. Only a simultaneous redefinition 
of the Breit-Wigner parameters for the higher-mass scalar Kq{14:30) and the 
inclusion of k gives an acceptable fit with a x^/d.o.f of 1.17, corresponding to a 
6.8% C.L. The mass and width of this effective K^{U30) are 1461 ± 4 MeV/c^ 
and 177 ± 8 MeV/c^, respectively; they can be compared to the PDG values 
of mass 1412 ± 6 MeV/c^ and width 294 ± 23 MeV/c^, and to the K-matrix 
fit values of mass 1408 MeV/c^ and width 220 MeV/c^. The parameters for 
the K in a fit with a free -^"0(1430) are mass 856 ± 17 MeV/c^ and width 
464 ± 28 MeV/c^. In Table El fit fraction^ phases and coefficients of the 
various amplitudes from the fit are reported, along with masses and widths 
for K and ii^^(1430). Coefficients refer to Bose-symmetrized, normalized am- 
plitudes, both for resonant and non-resonant states. Masses and widths can 
be compared with previous determinations from E791 [19] and BES [20]. All 



The quoted fit fractions are defined as tfie ratio between ttie intensity for a singie 
ampiitude integrated over tfie Dalitz plot and that of the total amplitude with all 
the modes and interferences present. 
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channel 


fvaction (%) 


phase 6i (dcg) 


coefficicTit 


non — resonant 


zy. r ± 4.5 


00 ir 1 A 

oZo ± 4 


1.4/ ± U.ii 




± i.o ± z.i (^see text j 


10110 
zt Z zt i.Z 


± U.Uo ± U.Uo 




io. / ± u.y 


U (nxeaj 


l(nxeaj 




± U.b ± U.O 






K (i4iUj7r^ 


no 1 n 1 
U.Z ± U.i 


350 ± 34 


n 1 1 n no 

U.lz ± U.U3 




± U.i ± U.U4 


1 1 ^ 1 1 r 
± i / ± i5 


± U.UU3 ± U.Ui 




1 Q _i_ n /I 
i.o ± U.4 


3 ± ( 


U.3d ± U.U4 




1 n 1 n 
± U.2 ± U.O 


± 4 ± 


1 n no -1- n nQ 
± U.Uz ± U.Uo 


A2 (i4oU jvr 


U.4 ± U.UO 


oiy ± 


n 1 7 ^ n ni 
U.i / ± U.Ui 




± U.U4 ± U.Uo 


1010 

± Z zh z 


1 n ni ^ n ni 
± U.Ui ± U.Ui 




i ( .0 It i.o 


OD It 


1 1 Q ^ n nc: 
i.io zt U.UO 




± 0.8 ± 0.4 


± 2 ± 1.2 


± 0.01 ± 0.02 




22.4 ± 3.7 


199 ± 6 


1 28 ± 10 




± 1.2 ± 1.5 (see text) 


± 1 ± 5 


± 0.015 ±0.04 




mass (MeV/c^) 


width (MeV/c2) 




Kg* (1430) 


1461±4±2±0.5 


177±8±3± 1.5 






856 ±17±5±12 


464 ± 28 ± 6 ± 21 





Table 2 

Fit fractions, phases, and coefficients from the isobar fit to the FOCUS — > 
-fC^TT^vr^ data. The first error is statistic, the second error is systematic from the 
experiment, and the third error is systematic induced by model input parameters 
for higher resonances. 



the results are consistent within the errors. Two systematic uncertainties are 
quoted in Table [2]: the first includes our split sample and fit variant compo- 
nents of Section 3.3, the second reflects uncertainties due to higher resonance 
modeling. In particular, this analysis is performed with a fixed radius of 5 
GeV~^ and 1.5 GeV~^ in the Blatt-Weisskopf D-meson and resonance form 
factors, respectively [19]. To cover the large range of uncertainties we found 
in literature on these parameters, we estimated the systematic effect by vary- 
ing them. Furthermore, for J > resonances we assumed central values for 
masses and widths from the PDG [18]. For those parameters which are not ac- 
curately determined we allowed a ±1 a variation and estimated the systematic 
uncertainty by the corresponding spread in the results. An analogous study is 
performed for the K-matrix fit and the corresponding uncertainties are quoted 
in Table [3] and HI Systematics for the k and non-resonant components have 
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Figure 3. Top: Dalitz plot projections with our isobar fit superimposed. The back- 
ground shape under the signal is also shown. Bottom: The adaptive binning scheme. 



a very large dependence on whether or not Gaussian form factors, suggested 
in [21], are included in the fit. Our systematic uncertainties do not account for 
the effect introduced by adding form factor for the scalar resonances. There 
is still some controversy on whether form factors are even needed for these 
scalars [22]. We find that adding in the Gaussian form factors does not im- 
prove our fit quality (3% C.L versus 6.8% C.L.). In this case, our results be- 
come consistent, within the errors, with Model C in [19], and changes the k and 
non-resonant fit fractions to (40.7 ± 6.2) % and (17.5 ± 4.1 )% respectively. 
Dalitz plot projections with fit results and the corresponding adaptive-binning 
scheme are shown in Fig. [31 
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3.5 K-matrix model results 



In the K-matrix fit to the K^n^n^ decay the free parameters are 

amplitudes and phases (oj and 5i) for vectors and tensors, and the P-vector 
parameters for scalar contributions. Kn scattering determines the parameters 
of the K-matrix elements and these are fixed inputs to this D decay analysis. 
Table [3] reports our K-matrix fit results. It shows quadratic terms in (-Pi/2)i 
are significant in fitting data, while in both (^1/2)2 and P3/2 constants are 
sufficient. 





coefficient 


phase (deg) 




= 3.389 ± 0.152 ± 0.002 ± 0.068 


e = 286 ±4± 0.3 ±3.0 


ClO 


= 1.655 ± 0.156 ± 0.010 ± 0.101 


71 = 304±6±0.4±5.8 


Cll 


= 0.780 ± 0.096 ± 0.003 ± 0.090 




C12 = 


= -0.954 ± 0.058 ± 0.0015 ± 0.025 




C20 


= 17.182 ± 1.036 ± 0.023 ± 0.362 


72 = 126 ±3 ±0.1 ± 1.2 


C30 


= 0.734 ± 0.080 ± 0.005 ± 0.030 


73 = 211 ±10 ±0.7 ±7.8 


Total S-wave fit fraction = 83.23 ± 1.50 ± 0.04 ± 0.07 % 
Isospin 1/2 fraction = 207.25 ± 25.45 ± 1.81 ± 12.23 % 
Isospin 3/2 fraction = 40.50 ± 9.63 ± 0.55 ± 3.15 % 



Table 3 

S'-wave parameters from the K-matrix fit to the FOCUS — > K~tt~^'k'^ data. 
The first error is statistic, the second error is systematic from the experiment, and 
the third is systematic induced by model input parameters for higher resonances. 
Coefficients are for the unnormalized S-wave. 

The J > states required by the fit are listed in Table HI 

The S'-wave component accounts for the dominant portion of the decay (83.23± 
1.50)%, and is the sum of the 7 = 1/2 and 1 = 3/2 components, which, sepa- 
rately, account for the (207 ± 24)% and (40 ± 9)% of the decay, with -164 % 
from their interference. The large amount of interference between the 1 = 3/2 
component and the 7 = 1/2 component underscores its importance in our 
K-matrix fit. Because there are no 7 = 3/2 resonances, the 7 = 3/2 S-wave 
component has at most a slowly varying phase and amplitude and since we 
find little variation in the 7^3/2-vector phase as well, the 7 = 3/2 piece essen- 
tially plays a comparable role to the ~ 30% non-resonant component present 
in the isobar fit summarized by Table 2. A significant fraction, 13.61 ± 0.98%, 
comes, as expected, from 7^*(892); smaller contributions come from two vec- 
tors 7r*(1410) and 7r*(1680) and from the tensor 7^2(1430). It is conventional 
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± 0.004 ±0.05 


± 0.3± 1 


± 


0.010 ±0.012 


i^*(1410)7r+ 


0.48 ±0.21 


293 ± 17 





.188 ±0.041 




± 0.012 ±0.17 


± 0.4± 7 


± 


0.002 ± 0.030 



Table 4 

Fit fractions, phases, and coefficients for the J > components from the K-matrix 
fit to the FOCUS — > K^tt^tt^ data. The first error is statistic, the second error 
is systematic from the experiment, and the third error is systematic induced by 
model input parameters for higher resonances. 

to quote fit fractions for each component and this is what we do. Care should 
be taken in interpreting some of these since strong interference can occur. This 
is particularly apparent between contributions in the same-spin partial wave. 
While the total S'-wave fraction is a sensitive measure of its contribution to 
the Dalitz plot, the separate fit fractions for 1 = 1/2 and 1 = 3/2 must be 
treated with care. The broad J = 1/2 5- wave component inevitably interferes 
strongly with the slowly varying J = 3/2 S'-wave, as seen for instance in [23]. 
Fit results on the projections and the adaptive binning scheme are shown in 

Fig. a 

The fit x^/d-o.f is 1.27 corresponding to a confidence level of 1.2%. If the 
I = 3/2 component is removed from the fit, the x^/d.o.f worsens to 1.54, 
corresponding to a confidence level of 10~^. 

3. 6 Comparison and discussion of the results 

The isobar fit represents a good effective description of the data, as testified by 
the confidence level. However, simple Breit-Wigner forms have been used for 
both the broad scalars k and i^'o(1430), each with free mass and width with 
no reference to how these states appear in other Kn interactions. Elastic scat- 
tering provides independent information about these states, their shape and 
parameters and how they overlap. In particular, scattering data from LASS 
show that the phase of the S'-wave rises by no more than 100° from 825 to 
1450 MeV, while the isobar fit with its simple Breit-Wigner forms requires 
~ 180° change in the resonant 7 = 1/2 S'-wave. In Fig. E^) and b) we compare 
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Figure 4. Top: The Dalitz plot projections with the K-matrix fit superimposed. The 
background shape under the signal is also shown. Bottom: The adaptive binning 
scheme. 

the modulus and phase of the total S'-wave components from the isobar and 
K-matrix fits. They essentially agree, as expected, since they fit the same data. 
However, the physics description differs considerably in the J = 1/2 compo- 
nents. In Fig.Oc) and d) modulus and phase of the J = 1/2 S'-wave component 
from our K-matrix fit are shown. The motivation for the K-matrix fit was to 
bring consistency between the description of scattering and D-decay data. In 
such a formalism the poles of the S-matrix are process independent, and on 
the real energy axis the overlap of broad resonances is correctly described. 
The results of the K-matrix fit showed that such a consistent representation is 
equally possible, the global fit quality being indeed good. However, it deterio- 
rates at higher Kn mass. This is not surprising since our K-matrix treatment 
only includes two channels Ktt and Krj' . While we have reliable information 
on the former channel, we have relatively poor constraints on the latter. This 
means that as we consider Kn masses far above Kri' threshold, these inade- 
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Figure 5. Modulus and phase of the S-wave components resulting from the K-matrix 
and isobar fit. Bands represent the ±1 a statistical variation of the K-matrix results 
for the total component [ a) and b) plots], the 1 = 1/2 component [ c) and d) plots)] 
and the the / = 3/2 component [ e) and f)] plots). Continuous lines in a) and b) 
represent the central value of the isobar fit. Vertical dashed line shows the location 
of the Kr]' threshold. 



quacies in the description of the Krj' channel become increasingly important. 
This is expected to become worse as yet further inelastic channels open up. 
Consequently, improvements could be made by using a number of D-decay 
chains with Kir final state interactions and inputting all these in one com- 
bined analysis in which several inelastic channels are included in the K-matrix 
formalism. In the present single K^ti^ti^ channel, adding further in- 

elastic modes would be just adding free unconstrained parameters for which 
there is little justification. It is interesting to note that the adaptive binning 
scheme shows that both the K-matrix and the isobar fit are not able to repro- 
duce data well in the region at 2 GeV^, in the vicinity of the Kr]' threshold. It 
is also the energy domain where higher spin states live. Vector and tensor fit 
parameters in the two models are in very good agreement: we do not exclude 
the possibility that a better treatment of these amplitudes could improve the 
X^- Some isolated spots of high could be caused by an imperfect modeling 
of the efficiency as they are in the same regions in both fits. The FOCUS 
experiment has also studied the K^n^ amplitudes in the K^K^n^ 
decay. We note that the moduli of Fig. [5^) are very different than that pre- 
sented in Fig. 3a) of [24]. This underscores the fact that different decay modes 
of the same charm state can often have quite different shapes for the moduli 
of their amplitudes, reflecting the differing production dynamics encoded in 
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the P-vectors. 



4 Conclusions 

The analysis of the — > K^ti^ti^ decay with the two models described in 
this paper reveal quite different features. The isobar model with its Breit- 
Wigner representation for all states requires both a k and a ii'Q(1430) whose 
parameters are not what elastic scattering would require. However, as already 
indicated, such Breit-Wigner parameters are effective rather than genuine 
pole positions. In contrast, the K-matrix fit has built in consistency with Ktt 
scattering. Moreover, this agrees with the results of our FOCUS experiment in 
semi-leptonic Z)-decay [25-27], where we have shown that the S'-wave phase 
agrees with that of the LASS experiment. The hypothesis of the two-body 
dominance, which has already been tested in other charm meson decays, is 
consistent with our results for the high-statistics K^-k^-k^ . Previous 

analyses [28] did not allow for the two different isospin states, and compared 
the 7 = 1/2 scattering phase with the global scalar phase. Comparison with 
our global F-vector phase is shown in the left plot of Fig. [^p] . A feature of the 
K-matrix amplitude analysis is that it allows an indirect phase measurement of 
the separate isospin components: it is this phase variation with isospin 1 = 1/2 
which should be compared with the same 1 = 1/2 LASS phase, extrapolated 
from 825 GeV down to threshold according to Chiral Perturbation Theory. 
This is done in the right plot of Fig.[6l As explained in Section 3.2, in this model 
[7] the P-vector allows for a phase variation accounting for the interaction with 
the third particle in the process of resonance formation. It so happens that 
the Dalitz fit gives a nearly constant production phase. The two phases in 
Fig. [UJo) have the same behaviour up to ~ 1.1 GeV. However, approaching 
Kfj' threshold, effects of inelasticity and differing final state interactions start 
to appear. The difference between the phases in Fig. [6^) is due to the 7 = 3/2 
component. 

These results are consistent with Kn scattering data, and consequently with 
Watson's theorem predictions for two-body Ktt interactions in the low Kir 
mass region, up to ~ 1.1 GeV, where elastic processes dominate. This means 
that possible three-body interaction effects, not accounted for in the K-matrix 
parametrization, play a marginal role. 

Our results for the total S'-wave are in general agreement with those from 
the E791 analysis, in which the S-wave modulus and phase were determined 

^ Phases determined from scattering are absolute. Those from the K^tt'^tt'^ 
Dalitz analysis are relative. We are free to raise the phases of Fig. to be zero at 
threshold. 
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1.25 1.5 

(a) M(K7i) (GeV/c') 



.25 1.5 

(b; M(Kn) (GeV/c^ 



Figure 6. Comparison between the LASS 1 = 1/2 phase + ChPT (continous hne) 
and the F-vector phases (with ztlcr statical error bars); a) total F-vector phase; b) 
7 = 1/2 F-vector phase. Vertical dashed line shows the location of the Krj' . 



in each Kit slice [28], [29]. Edera and Pennington [23] were able to separate 
this total 5- wave into 7 = 1/2 and 7 = 3/2 components, making the strong 
assumption that the equivalent of the P-vector has a constant phase. Here 
we have relaxed this assumption, and find a slightly different separation of 
7 = 1/2 and 3/2 components, but the trend with energy and relative phases 
are broadly consistent. 

What does this analysis contribute to the discussion of the existence and pa- 
rameters of the k7 We know from analysis [30] of the LASS data (which in 
7r~7r+ scattering only start at 825 MeV) there is no pole, the k(900), in its 
energy range. However, below 800 MeV, deep in the complex plane, there is 
very likely such a state. Its precise location requires a more sophisticated ana- 
lytic continuation onto the unphysical sheet than the K-matrix representation 
provided here. This is because of the need to approach close to the crossed 
channel cut, which is not correctly represented for a robust analytic contin- 
uation. However, our K-matrix representation fits along the real energy axis 
inputs on scattering data and Chiral Perturbation Theory in close agreement 
with those used in the analysis by Descotes-Genon and Moussallam [31] that 
locates the k with a mass of (658 ± 13) MeV and a width of (557 ± 24) MeV 
by careful continuation. These pole parameters are quite different from those 
implied by the simple isobar fits presented here, and by E791 in [19]. What 
we have shown is that whatever k is revealed by our — >■ K'h^tt^ results, 
it is the same as that found in scattering data. Consequently, our analysis 
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supports the conclusions of [31] and [32]. 

We have seen that i^-decay can teach us about Kn interaction much closer 
to threshold than the older scattering results. This serves as a valuable check 
from experiment [33] of the inputs to the analyses of [31] and [32] based largely 
on theoretical considerations. Indeed, more complete experimental insight into 
the Kn interaction will be provided by the full range of hadronic and scmi- 
leptonic /^-decays to come from 5-factories. Our results show that the dy- 
namics of the ir~7r^7r+ final state is dominated by two-body Ktt interactions 
up to 1.1 GeV as determined by scattering experiments. 
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